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Abstract 



We examine a model of non-self-avoiding, fluctuating surfaces as 
a candidate continuum string theory of surfaces in three dimensions. 
This model describes Dynamically Triangulated Random Surfaces em- 
bedded in three dimensions with an extrinsic curvature dependent 
action. We analyze, using Monte Carlo simulations, the dramatic 
crossover behaviour of several observables which characterize the ge- 
ometry of these surfaces. We then critically discuss whether our ob- 
servations are indicative of a phase transition. 
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1 Introduction 



In this work, we investigate a theory of fluid, fluctuating random surfaces em- 
bedded in three dimensions. Various theories of fluctuating surfaces (string 
theories) arise in the description of many physical systems. Among these 
are the 3-dimensional Ising model, the strong interaction (QCD) in the in- 
frared limit and fundamental unifled theories of all interactions including 
gravity. Natural biological membranes, such as lipid bilayers, and artiflcial 
membranes, such as micelles and vesicles, also form a rich class of fluctuating 
surfaces together with interfaces, such as those in microemulsions, between 
two distinct three-dimensional bulk phases [1]. These latter systems are fluid 
because their component 'molecules' are loosely bound. Their constituents 
are arranged so that their net surface tension nearly vanishes; thus these 
membranes undergo large thermal fluctuations. These biological and chemi- 
cal membranes exhibit self-avoidance, which we do not take into account in 
our simulations. 

Just as fleld theories are described by sums over paths, string theories are 
formally characterized by a functional integral Z which can be written as a 
sum over surfaces weighted by exp(— 5"). Here S denotes the action associated 
with a particular surface. Note that we work in Euclidean space. To write 
down our action we introduce an explicit parametrization of a generic surface 
M. in i?^ with coordinates {(J\^(J2) and the embedding X^{ai). \x runs from 
1 to 3 (since we only study the case of a 3d embedding space) . The induced 
metric (the puUback of the Euclidean metric via the embedding) is given 



We will use Greek letters for the embedding space indices; Roman letters 
label the coordinate basis on each surface. A particularly natural choice for 
the action is the Nambu-Goto action, which is proportional to the area of 
the surface in the embedding space, and is given by 



This action is difficult to quantize, though, since it is non-polynomial. A 
further complication is that the measure in the path-integral constructed 
with this action is quite subtle, especially when the theory is discretized. By 



by 



(1) 




(2) 
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introducing additional degrees of freedom, in the form of an intrinsic metric 
Qij on each surface, one can write down the Polyakov action ^. 



This action has proven to be much more tractable analytically. Note that 
it possesses a large degree of gauge freedom, associated with its invariance 
under reparametrizations of the metric and the intrinsic surface coordinates. 
Often, one gauge-fixes the metric to the form gij = exp(p)5.y; p is referred 
to as the Liouville mode. The Polyakov action is in fact independent of p ^ , 
but this 'Weyl' symmetry is anomalous in general. The anomaly is expressed 
in terms of the central charge c. By carefully gauge-fixing the entire path- 
integral one can derive the anomaly condition 



The first term depends on the space in which the surfaces are embedded, or, in 
other words, on the particular field theory that lives on the two-dimensional 
surface. For instance, the central charge associated with D bosons is D (1 per 
boson); this characterizes a string embedded in ^ The —26 arises from the 
anomaly contribution from ghosts needed for gauge-fixing reparametrization 
invariance. Note that when D = 26, Cp, the Liouville central charge, vanishes 
and the Liouville mode effectively decouples from the theory. For D < 26 
one can integrate the matter (X^) fields out and obtain an effective Liouville 
theory. The Liouville theory is highly non-linear; furthermore, its analysis 
is complicated by the fact that the two-dimensional metric on the surface in 
which the Liouville field lives depends on the field itself! Nevertheless, with 
a few general assumptions about the path integral measure, the spectrum 

"'^ An additional term, which is a topological invariant proportional to the Eulcr character 
of the surface is added to both the Polyakov and Nambu-Goto actions. The coefficient of 
this term is referred to as the 'string coupling'. It is the only dimensionless parameter in 
the theory; therefore it serves as the pcrturbativc expansion parameter for string theory. 

^It follows that substituting the solution to the equations of motion back into the 
Polyakov action yields the Nambu-Goto action. At the quantum level, though, it has only 
been shown that these actions are equivalent for D = 26 [2]. In lower dimensions this 
equivalence has been questioned, for example, by Distler [3]. The work of Polchinski and 
Strominger also suggests that there are alternate quantizations [4]. 

^ A fermion has a central charge of 1/2, so that by studying an Ising model on random 
fluctuating surfaces, we are effectively studying strings imbedded in 'fractional' dimensions. 




(3) 



Cmatter + Cp — 26 — 0. 



(4) 
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of the theory at each order in perturbation theory can be computed when 
Cmatter < 1- Additional exponents, such as the Hausdorff dimension, which 
characterize both the intrinsic and extrinsic geometry of these surfaces, have 
also been computed from Liouville theory. These theories are often analyzed 
on surfaces of fixed area; in this case, the string susceptibility is determined 
by the dependence of the generating functional for surfaces of genus h on the 
area A 

\n{Zh) = exp(/Xe^)^^''-^ (5) 

where Hc is a non-universal constant. An analysis of the Liouville model 
yields 

7fe = 2 - ^ (25 - Cmatter + -^(1 - Cinatter)(25 - Cmatter))- (6) 

These computations break down in the regime 1 < Cmatter < 25, where expo- 
nents such as the string susceptibility given in (6) become imaginary. 

Further progress in understanding string theory in a low number of di- 
mensions has been made by mapping the sum over surfaces onto an integral 
over matrices. This mapping is realized by replacing continuum surfaces by 
their discrete cellular decompositions, such as their triangulations. In this 
construction, each triangle is equilateral with area 1 in the intrinsic metric; 
the coordination number at each vertex determines the intrinsic curvature 
of the surface. The coordinates i label the vertices of the triangulation. 
Then the discrete analogue of the intrinsic metric is the adjacency matrix 
Cij whose elements equal 1 if i and j label neighbouring nodes of the trian- 
gulation, and vanish otherwise. Two-dimensional diffeomorphism invariance 
reduces, at this discrete level, to the permutation symmetry of the adjacency 
matrix. One of the keys, in fact, to the power of this construction is the 
preservation of this symmetry. The triangulation of a surface of genus h is 
then dual to a phi-cubed diagram of genus h. The large N ('t Hooft) expan- 
sion of an integral over N x N Hermitian matrices generates these Feynman 
graphs. In this case the area of the surfaces is not fixed, since all graphs 
of fixed genus are summed over. The action then acquires a contribution 
(from Legendre transformation) proportional to the product of the cosmo- 
logical constant with surface area A. Orthogonal polynomial techniques can 
then be applied to analyze the matrix integral. Indeed the theory can be 
exactly solved in the 'double scaling limit' [5] in which the string coupling. 
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cosmological constant and matter couplings (which are determined by the 
matrix 'potential') are tuned together. These solutions are exact, and thus 
include non-perturbative information about these string theories, albeit with 
ambiguities in certain cases. Unfortunately, the matrix models that represent 
theories with Cmattcr > 1 ai'e too difficult to solve exactly. 

It has generally been suspected that these analytic techniques fail for 
Cmatter > 1 because the string theory becomes pathological. In this regime, 
a negative mass-squared particle, which is sometimes referred to as the 
'tachyon', comes on shell ^. It is thought that its presence might engender the 
proliferation of long tubes with thickness of the scale of the ultraviolet cutoff 
on surfaces that dominate the string functional integral. Gates, by analyzing 
the Liouville action, has put forth another perhaps related explanation of this 
pathology [7]. He points out that in the presence of vortex configurations 
of the form p ~ — ln(r), reducing the world-sheet cutoff to zero no longer 
causes the spacetime cutoff to vanish. He then shows that these vortices have 
negative free-energy when Cmattcr > 1 and proliferate to the extent that the 
Liouville partition function becomes ill-defined with a finite spacetime cutoff 
and zero cosmological constant. For positive cosmological constant, it then 
appears plausible that these vortex configurations would still appear, though 
perhaps the partition function will be finite. The predominant surfaces in 
the functional integral would thus be subject to large fiuctuations of their 
internal geometry. Recent simulations of multiple Potts models with Cmatter 
less than, equal to, and greater than 1 coupled to gravity have been per- 
formed [8]. These simulations, however, do not show any dramatic changes 
in the behaviour of the internal geometry of the dominant surfaces as Cj„attcr 
is increased beyond 1. More analytic and numerical work clearly needs to be 
done to determine what exactly happens as Cmatter becomes greater than 1. 

Monte Carlo simulations of strings embedded in flat space for D > 1 do 
indicate that these theories fail to describe the fluctuations of two dimensional 
smoothly embedded surfaces in the continuum limit. The normals to the 
surfaces dominating the simulations are uncorrclatcd over the distance of 
a few lattice spacings. These surfaces also appear to have a large (greater 
than 8), or perhaps infinite, extrinsic Hausdorff dimension; they resemble 

*More precisely, one encounters these instabilities in Liouville theory when the quantity 
c — 24 A > 1, where c denotes the central charge of the matter theory which describes the 
embedding of the surfaces and A is the conformal weight of the lowest weight state in this 
theory [6]. Since here we are considering flat space, c = D and A = 0. 
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branched polymers ^ . Such configurations should not describe, for instance, 
the domains of the 3d Ising model or QCD strings. We would thus like to 
find a modified string theory that is dominated by smoother surfaces. 

The tachyon, and apparently these related instabilities, can be eliminated 
in particular cases by introducing fermionic coordinates and supersymmetry 
on the worldshcct, and implementing an appropriate projection of states. 
The fermions, presumably, effectively smooth out the surfaces. This would 
be consistent with what has been observed for one-dimensional geometries; 
the random walk of a spin one-half particle has Hausdorff dimension one and 
thus appears to be smooth [9]. Many authors have proposed an alternative 
modification of the string action [10, 11, 12, 13] via the addition of a term 
that directly suppresses extrinsic curvature^. We shall examine this class of 
theories in this talk. 

To characterize the geometry of our surfaces further, we associate with 
each point in our generic surface M. tangent vectors (tf^TAi) and a normal 
vector n^^ETM.^. The extrinsic curvature matrix Kij (the second fundamen- 
tal form) can be defined by 

diU^ = -Kijt^^ . (7) 

The eigenvalues of this matrix are the inverses of the radii of curvature of M.. 
One usually describes the geometry of these surfaces in terms of the mean 
curvature [15, 16] 

H = h^^K.j , (8) 

and the Gaussian curvature 

K = e'h^'K.^Kki . (9) 

One can show that the Gaussian curvature can be computed solely from 
the metric /ijj, while the mean curvature depends explicitly on the embedding 

^As above, singular configurations also dominate the Gaussian theory, which is essen- 
tially a theory of free random walks, rather than surfaces. 

^In fact, integrating the fermions out of the Green-Schwarz superstring yields an action 
similar to the one we consider, but with the addition of a complex Wess-Zumino type term 
[14]. 
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Our lattice model is constructed by triangulating each surface, as we dis- 
cussed above in the context of matrix models. Each node of the triangulation 
is embedded in by the functions ; i labels the ith node and /j, runs from 
1 to 3. We also associate a normal vector {n^)f, with each triangle (indices 
with hats label the triangles). We shall study the theory defined by the 
action 



S^Sg + XSe = Y: CijiXt - X^f ^XY^C^'il-nl nf ). (10) 



This model has been examined in references [17, 18, 19, 20, 21, 22, 23]. For 
A = this is simply a discretization of the Polyakov string action. The 
final term, which depends on the discretized extrinsic curvature, introduces 
a ferromagnetic interaction between surface normals, which one might hope 
would cause smoother surfaces to dominate the partition function. We would 
like to know if there is a smooth phase and a phase transition (at some finite 
Ac) between this phase and the crumpled phase observed at A = 0. If this 
were so, an interesting continuum limit of this lattice model could perhaps be 
constructed at this phase transition point, yielding a new continuum string 
theory. 

The action we simulate is in fact a particularly natural discretization of 



Note that the second term in the action is manifestly positive, Weyl and 
reparamctrization invariant, and that A is a dimcnsionless coupling. So, it 
is not clear whether or not this term is relevant. Mean field (large D) and 
perturbative RG calculations have been performed using similar actions, such 
as 



in the regime in which the string tension //q is small (unlike the usual particle 
physics limit of string theory, which is characterized by large //q)- After 
integrating out fluctuations of the embedding X'^ between momentum scales 
A and A, it is found that the renormalization of the extrinsic curvature 
coupling is given to one-loop order by 



k,i,fj, 




(11) 




(12) 



m \ I, da 3 2 



(13) 
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so that a is driven to infinity in the infra-red. This theory thus exhibits 
asymptotic freedom. Surfaces are smooth (the normals are correlated) below 
a persistence length [24] 




(14) 
can be 





gained by observing that this theory is similar to the 0(3) sigma model, 
which is asymptotically free [9]. The normals to M arc the analogues of 
0(3) vectors, though in this case they are constrained to be normal to a 
surface governed by the action (11). 

The analytic results, therefore, do not indicate that we should anticipate a 
finite A phase transition. Note, though, that the RG calculations are based on 
a Nambu-Goto type action although we simulate an extension of the Polyakov 
action. Since the two actions are not clearly equivalent, 'particularly when 
extrinsic curvature dependent terms are added, we cannot simply assume 
that these analytic and our numerical results should agree. 



We have considered triangulations with the topology of the torus, to minimize 
finite size effects. The above action was used, with the BRST invariant 
measure utilized also by Baillie, Johnston and Williams [18], so that 



where D — 3, Qi is the connectivity of the ith vertex and T(l) refers to the 
set of triangulations of genus 1. We used the standard Metropolis algorithm 
to update our configurations. To sweep through the space of triangulations 
we performed flips on randomly chosen links. Flips were automatically re- 
jected if they yielded a degenerate triangulation. After a set of 3M flips 
was performed, 3M randomly selected embedding coordinates were updated 
via random shifts from a fiat distribution. Most of the Monte Garlo simula- 
tions were performed on HP-9000 (720 and 750 series) workstations; we also 
collected some data by simulating lattices on each of the 32 nodes of a GM-5. 



2 The Simulation 



(15) 



GeT(l) lJ,,i i 
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N 




Ac 


36 


3.484(8) 


1.425(35) 


72 


4.571(15) 


1.410(15) 


144 


5.37(14) 


1.395(30) 


288 


5.55(7) 


1.410(25) 


576 


5.81(17) 


1.425(30) 



Tabic 1: The maximum of the specific heat and its position, with errors, for 
different lattice sizes. 

We ran on lattices ranging in size from A'^ = 36 to 576 {N signifies the 
number of vertices) with 4 to 7 different values of A for each N. Most of 
the data was this data was taken in the region A e (1.325, 1.475). For small 
A^, the runs consisted of 3 x 10^ sweeps, while we performed longer runs (of 
up to 27 X 10^ sweeps for = 576) for larger lattices, because the auto- 
correlation times for our simulations were very large. (The correlation time 
for the radius of gyration was greater than 10^ sweeps for N — 576!) To 
determine our observables as a function of A we used a histogram reconstruc- 
tion procedure. We patched different histograms by weighting them with the 
associated statistical indetermination (which was estimated by a jack-knife 
binned procedure). Various consistency checks indicate that this procedure 
is very reliable. 

3 Observables 

We measured the edge action Se and the associated specific heat C(A) = 
jj-{< S"^ > — < Se >^)- In Fig. 1 we plot the specific heat curve (con- 
structed via the histogram procedure) and we tabulate its maximum and 
peak position for various lattice sizes in Table 1. 

We see that the specific heat peak grows vigorously with N for small 
lattices, but that this growth quickly levels off for larger N. These observa- 
tions agree fairly well with previous work [20, 21, 22]. For the larger lattices 
it appears that the peak position shifts very slowly towards higher values 
of A, though this increase is not statistically significant. The shape of the 
peak does not change dramatically with N; it narrows perhaps a bit between 
N = 144 and N = 576. 
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To determine how the mean size of the dominant surfaces depend on A, 
we measured the squared radius of gyration Rq] 

Rg ^ ^T^i^t - xu\ (16) 

where the com subscript refers to the center of mass of the surface. By 
measuring the N dependence of the gyration radius, one can extract a value 
for the extrinsic Hausdorff dimension, which is given by 

Rq N" . (17) 

Wc plot Rq in Fig. 2(a), clearly the size of the mean surface size increases 
dramatically for A near 1.4. In Fig. 2(b), we plot the effective Hausdorff 
dimension, given by 

In the large A limit i/ — > 1 and dgxtr — > 2, as expected for fiat surfaces. 
In the low A limit dextr becomes very large. In the pseudo-critical region z/ is 
a linear function of A. Curiously enough, the latter curve yields a Hausdorff 
dimension of 4, a value characteristic of branched polymers, near the location 
of the specific heat peak. This value is not particularly reliable though be- 
cause of finite-size effects and also because it changes rapidly in this region. 
In ref. [22] a value compatible with ours {Dh{Xc) > 3.4) is quoted for the 
critical theory. We stress however (and also here we are in complete agree- 
ment with [22]) that the Hausdorff dimension in the pseudo-critical region 
depends heavily and quite unusually on N. 

In both the high and low A regions finite size effects are quite small 
(compatible with zero to one standard deviation). In the pseudo-critical 
region, on the contrary, finite size effects are large. This effect cannot be 
explained by the shift in A which one gets from the shift of the peak of the 
specific heat, which is far too small. 

We also measured the magnitude of the extrinsic Gaussian curvature. 



/ I ^ I \/r^> given by 
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Here (pj denotes the angle subtended by the jth triangle at the ith vertex. 
This quantity, plotted in Fig. 3, measures the magnitude of the deficit angle 
in the embedding space averaged over all vertices. Note that the mean Gaus- 
sian curvature decreases rapidly in the neighborhood of A = 1.4, indicating 
that a sharp crossover is occuring in this system. From this plot we can see 
that finite size effects increase with A. They do not appear to peak in the 
region about A = 1.4 as one might expect for a typical phase transition. 
The magnitude of the intrinsic Gaussian curvature, | 7^ | , given by 

l^h^E|6-gd , (20) 

i 

is shown in Fig. 4. When the intrinsic and extrinsic metrics are equal, 
the intrinsic and extrinsic deficit angles are identical and K = R/2. Both 
extrinsic and extrinsic curvatures behave in a qualitatively similar manner; 
I TZ I drops off rapidly, just as | /C | does. Through the peak region, though, 
I /C I decreases by about a factor of 5 while | TZ \ diminishes to only about 
.6 of its value on the left-hand side of the peak. Since the action explicitly 
suppresses mean curvature, and the mean and extrinsic Gaussian curvature 
are closely related (for instance, H"^ > y), we would expect that for large 
A extrinsic fluctuations would be suppressed much more than fluctuations of 
intrinsic geometry. 

The question of the equivalence of the Nambu-Goto and Polyakov actions 
motivated us to study the correlations between intrinsic geometry (which is 
not introduced independently in the Nambu-Goto formulation) and extrinsic 
geometry. We measured the quantity which we refer to as K. *TZ 

/^^ E.(2.-E,4)(6-..) 
'JSK'JR' yE<(27.-E5«yE<(6-%)^ 

This quantity is 1 when the metrics are equal, if they are un- correlated, 
and negative when these curvatures are anti-correlated. The plot of /C * 7?. in 
Fig. 5 indicates that intrinsic and extrinsic geometry are strongly correlated 
for small A, but as one passes through the peak region they become decor- 
related. This is not particularly surprising, given that the action directly 
suppresses only extrinsic fluctuations. Note that RG calculations based on 
the Nambu-Goto action plus an extrinsic curvature term perturb about a 
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background that is both intrinsically and extrinsically flat. Given the ob- 
served decorrelation between intrinsic and extrinsic geometry, we would not 
anticipate that this background appears in the low- temperature limit of the 
model which we simulate. 

We also measured various other observables which characterize both the 
intrinsic and extrinsic geometry of these surfaces. These measurements are 
discussed in another write-up of this work [23]. They all exhibit sharp 
crossover behaviour in the region near A — 1.4. We found that the auto- 
correlation times of these observables grew rapidly as A increased, but we 
did not note any maximum in these times in the region about A = 1.4. 

The crossover behaviour became also quite apparent when we examined 
typical snapshots of our simulated surfaces for various values of A. In Figs. 
6 (a)-(d) we present pictures of typical surfaces for A = .8, 1.3, 1.5, and 2.0. 
Note that the surfaces rapidly become smoother and the normal-normal cor- 
relation length increases significantly as one passes from the second to the 
third of these pictures, which correspond to only slightly different values of 
A. 

4 Interpretation 

This model of crumpled surfaces appears to exhibit sharp crossover behaviour 
in the region around A = 1.4. The sharp change in the magnitude of the 
Gaussian curvature, the radius of gyration and other observables indicates 
that the normals acquire long-range correlations, up to the size of the systems 
we examine. The zero string tension measurement of [22] also shows that the 
disordered regime differs from the regime in which the surfaces are ordered 
(up to scale of the lattices that are simulated) by only a small shift in A. This 
evidence might indicate the presence of a phase transition at this point. Since 
the peak growth rapidly diminishes for large N, such a phase transition would 
likely be higher than second order. Still, the apparent absence of diverging 
correlation times and, in some cases, increasing finite size effects in the peak 
region leads us to question whether we are actually observing a typical phase 
transition. 

There are indeed other possible interpretations of our data. Note that 
the surfaces which we simulate are quite small. For instance, if the surfaces 
in our simulations had an intrinsic dimension of 2.87 (characteristic oi D — 
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gravity), they would have roughly a linear size of fewer than 10 lattice 
spacings 

Perhaps the simplest alternative explanation for the presence of this peak 
is suggested by the arguments of KroU and Gompper [25]. They argue that 

the peak occurs when the persistence length of the system approaches the 
size of the lattice {^p ~ Nd); d denotes the intrinsic Hausdorff dimension. 
Fluctuations on a larger scale become more important. When this scale is 
greater than the lattice size these fluctuations are suppressed. Thus one 
might surmise that the specific heat will drop for large A. Typically, the 
persistence length grows as ~ exp(CA); C is inversely proportional to 
the leading coefficient of the beta function. We would expect that the peak 
position should shift to the right with increasing N in this scenario as 



Quite a large value of C is needed to explain the rapid crossover; roughly 
values of C ~ 10, dintr ~ 3 are more or less consistent with the magnitude 
of the peak shift and crossover width. The RG calculations using different 
forms of the action yield C = ^ (sec equation 14), but this may not apply 
to the action we simulate. This reasoning also indicates that the peak should 
widen as the lattice size increases; we do not observe this at all. It seems 
plausible though that these arguments, based only on the leading term of the 
high A expansion, are too naive. 

An alternative scenario, which builds on the ideas in the above paragraph, 
is suggested by the tantalizing similarities between the results of our fluid 
surface simulations and what has been observed for the d = 4 SU{2) Lattice 
Gauge Theory [26] and for the c? = 2 0(3) model. 

The 0(3) model, which is thought to be asymptotically free, exhibits a 
specific heat peak near /5 = 1.4 (first measured via Monte Carlo simulations 
by Colot [27]). The origin of this peak is understood [28]; it is due to the fluc- 
tuations of the sigma particle, a low-mass bound state of the massless 0(3) 
pions. The sigma induces short-range order and contributes to the specific 
heat as a degree of freedom only at high temperatures (when the correlation 
length in the system becomes smaller than its inverse mass). The peak thus 

''Of course, our lattices are too small, by one or two orders of magnitude, to really 
exhibit a convincing fractal structure. 





(22) 



dC 
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occurs at the beginning of the crossover regime, when the correlation length 
is several lattice spacings. 

According to the low temperature expansion, the correlation length grows 
as ,^ ~ exp(2nP) / p. Thus one would expect a fairly rapid crossover in the 
0(3) model; the correlation length should increase by roughly a factor of 9 
when /? is shifted by about 0.35. Such a crossover is indeed observed, though 
it is not so apparent that it is as dramatic as the crossover behaviour observed 
for fluid surfaces. ^ 

Recent simulations of the 0(3) model [29] indicate that the specific heat 
peak grows significantly when the lattice size L is increased from 5 to 15, 
and that virtually no growth in peak height is evident as L is increased 
further up to 100. Furthermore, the peak position shifts to the right as 
L grows and then appears to stabihze for large L. This is more or less 
what we observe in our simulations of fluid surfaces, on lattices of small size. 
We point out these similarities largely to emphasize that there does exist 
an asymptotically free theory (with low mass excitations) which exhibits 
crossover behaviour qualitatively similar to that observed in our simulations. 
The analogy is perhaps deeper, though, since the fluid surface action (with 
extrinsic curvature) in certain guises looks like a sigma model action. It 
would not therefore be so surprising from this point of view to find a sigma 
particle in these theories, perhaps associated with (n^ — 1), in which n denotes 
the unit normal to our surfaces. 

Another additional possibility is that fluctuations of the intrinsic geome- 
try (the Liouville mode) are responsible for short-range order and contribute 
to the specific heat peak. 

5 Conclusion 

We have introduced a model of fiuid random surfaces with an extrinsic cur- 
vature dependent action and explored its phase diagram. Unfortunately, we 
have been unable to determine if our model undergoes a phase (crumpling) 
transition at finite coupling. We have observed dramatic crossover behavior 
for particular observables in our Monte Carlo simulations, but on the other 

^To compare quantitatively the width of the crossover regimes for these two models it 
would be necessary to measure a correlation length (perhaps extracted from the normal- 
normal correlation function) in these random surface simulations. 
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hand, the correlation times and certain finite-size effects do not behave as 
one would expect in the presence of a phase transition. The behavior of 
other lattice models also indicates that it is possible that we are observing 
the effects of finite-mass excitations on small lattices, rather than a phase 
transition. We hope that future work will clarify this murky state of affairs, 
to determine if there indeed exists a crumpling transition for fluid surfaces. 
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7 Figure Captions 

Fig. 1 The specific heat C(A) as a function of A. As in all other pictures, filled 

circles and a dotted line correspond to iV = 144, crosses and a dashed 
line indicate N = 288, and empty squares and a solid line represent 
N = 576. 

Fig. 2a The gyration radius Rq defined in (16), plotted as in Fig. 1. 

Fig. 2b The effective inverse Hausdorff dimension u as a function of A, as de- 
fined in (17). The filled dots and the dashed curve are from a fit to 
the = 288 and N = 144 data, while the empty dots and solid curve 
represent the fit to A" = 576 and A^ = 288. 

Fig. 3 The extrinsic Gaussian curvature | /C | defined in (19), plotted as in 
Fig. 1. 

Fig. 4 The intrinsic curvature | TZ \ defined in (20), plotted as in Fig.l. 

Fig. 5 The intrinsic extrinsic curvature correlation, as defined in (21), plotted 
as in Fig.l. 

Fig. 6a A snapshot of a 576 node surface of toroidal topology for A = 0.8. 
Fig. 6b As in Fig. 6a for A = 1.3. 
Fig. 6c As in Fig. 6a for A = 1.5. 
Fig. 6d As in Fig. 6a for A = 2.0. 
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